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ABSTRACT 

We examine whether the spectral energy distribution of optical continuum emission of active galactic 
nuclei (AGNs) changes during flux variation, based on accurate and frequent monitoring observations 
of 11 nearby Seyfert galaxies and QSOs carried out in the B, V, and / bands for seven years by the 
MAGNUM telescope. The multi-epoch flux data in any two different bands obtained on the same 
night show a very tight linear flux to flux relationship for all target AGNs. The flux of the host galaxy 
within the photometric aperture is carefully estimated by surface brightness fitting to available high- 
resolution HST images and MAGNUM images. The flux of narrow emission lines in the photometric 
bands is also estimated from available spectroscopic data. We find that the non-variable component 
of the host galaxy plus narrow emission lines for all target AGNs is located on the fainter extension 
of the linear regression line of multi-epoch flux data in the flux to flux diagram. This result strongly 
indicates that the spectral shape of AGN continuum emission in the optical region 4400-7900A) 
does not systematically change during flux variation. The trend of spectral hardening that optical 
continuum emission becomes bluer as it becomes brighter, which has been reported by many studies, 
is therefore interpreted as the domination of the variable component of the nearly constant spectral 
shape of an AGN as it brightens over the non-variable component of the host galaxy plus narrow lines, 
which is usually redder than AGN continuum emission. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — galaxies: quasars: general 
— accretion, accretion disks — dust, extinction 



1. INTRODUCTION 

The UV-optical variability of type 1 AGNs such as 
Seyfert 1 galaxies and QSOs was found at the begin- 
ning of AGN studies in the 1960s and has been used as a 
powerful tool for understanding the nature of AGN cen- 
tral engines. The mechanism of their variability is under 
active discussion, and proposed models cover a diverse 
range of mechanisms including the instability of an accre- 
tion disk (Rees 1984; Kawaguchi et al. 1998), X-ray re- 
processing (Krolik et al. 1991), star collision (Courvoisier 
et al. 1996), multiple supernovae (Terlevich et al. 1992), 
and gravitational microlensing (Hawkins 1993). None of 
these models, however, successfully explained more than 
a few properties of UV-optical variability (Vanden Berk 
et al. 2004). 

In efforts to construct a successful model of variability, 
photometric and spectroscopic monitoring observations 
have shown that the amplitude of UV-optical variability 
correlates well with various quantities of AGN luminos- 
ity, central black hole mass, accretion rate, redshift, rest- 
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frame wavelength, and so on (Vanden Berk et al. 2004; 
Wold et al. 2007; Wilhite et al. 2008). Among these, a 
key relation that directly reflects the nature of the cen- 
tral engine is a blueing trend, in which AGNs become 
bluer when they brighten. Giveon et al. (1999) moni- 
tored 42 PG quasars in the BR bands for seven years 
and showed such a blueing trend for at least half the 
quasars in their sample. Webb & Malkan (2000) mon- 
itored 23 AGNs in their special blue, yellow, and red 
bands for three months and showed a larger amplitude of 
variability in the blue band than that in the red band on 
timescales of months, while they found no color depen- 
dence of the variability amplitude on shorter timescales. 
Vanden Berk et al. (2004) used the two-epoch multicolor 
data of about 25,000 quasars from the Sloan Digital Sky 
Survey (SDSS) and found an anti-correlation of variabil- 
ity amplitude with rest-frame wavelength especially at 
A < 4000A in the UV-optical region. 

The blueing trend, although observationally estab- 
lished, is known to admit a dual interpretation. Some 
authors claim a naive interpretation, so-called spectral 
hardening, such that the variable component becomes 
brighter and gets bluer (e.g., Wamsteker et al. 1990; 
Giveon et al. 1999; Webb & Malkan 2000; Vanden Berk 
et al. 2004; Wilhite et al. 2005), while others claim a 
different interpretation such that the variable component 
of constant blue color becomes brighter and increasingly 
dominates over the non- variable component of red color 
(e.g., Choloniewski 1981; Winkler et al. 1992; Winkler 
1997; Paltani & Walter 1996; Suganuma et al. 2006). 

The naive interpretation of spectral hardening was in- 
deed countered in a convincing way by Winkler (1997), 
who monitored 92 Seyfert 1 galaxies in the UBVRI 
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bands for three years and showed a nearly hnear flux 
to flux relationship in any two different bands obtained 
for each AGN on the same night. Using accurate and fre- 
quent BV monitoring data for a few Seyfert 1 galaxies by 
the MAGNUM telescope (Kobayashi et al. 1998a, 1998b; 
Yoshii 2002; Yoshii, Kobayashi, & Minezaki 2003), the 
linear relationship was undoubtedly confirmed by Sug- 
anuma et al. (2006) and Tomita et al. (2006). Winkler 
(1997) interpreted the linear relationship as indicating 
the constant color or constant spectral shape of the vari- 
able component, as originally proposed by Choloniewski 
(1981), with the assumption of non- variable component 
of the host galaxy located exactly on the fainter exten- 
sion of the straight line fitted to the flux data in the flux 
to flux diagram. 

The effect of host galaxy on the optical color vari- 
ability of AGNs has been discussed by many authors 
(Doroshenko & Lyutyi 1994; Giveon et' al. 1999; Webb 
& Malkan 2000; Trevese & Vagnetti 2002). Doroshenko 
& Lyutyi (1994) examined the correlation between the 
fluxes in two different bands of UBVR for 30 Seyfert 1 
galaxies and found, opposing to the Choloniewski's in- 
terpretation, too large a scatter of their data from the 
linear relationship to be consistent with the estimated 
errors for their data, and the host-galaxy component not 
on the fainter extension of the best-flt straight line for 
five out of their 30 targets. However, Winkler (1997) 
presented the much tighter linear relationship than the 
graphs displayed by Doroshenko & Lyutyi (1994), with 
the scatter consistent with the errors for most of their 
92 targets, hence they concluded that the large scatter 
found by Doroshenko & Lyutyi (1994) was likely caused 
by their underestimate of the measurement errors. 

Giveon et al. (1999) estimated the host-galaxy flux for 
ten targets in their sample from the HST WFPC2 images 
obtained by Bahcall et al. (1997), but did not explic- 
itly take into account its contribution in their discussion 
on the spectral hardening. Webb & Malkan (2000) esti- 
mated the host-galaxy flux for more than two thirds of 
their sample, but did not examine its effect for each of 
individual targets. In fact, they constructed a spectral 
model consisting of typical AGN and galaxy spectra, and 
examined the amplitude of multicolor variations. They 
assumed that the host-galaxy component contributes one 
third of the total flux in the yellow band allowing 0.5 
mag flux variation for the AGN component, as an ex- 
treme case for a large contribution of host galaxy, and 
found that the blue to red amplitude ratio was larger 
than that expected from the model for two out of their 
6 targets. Trevese & Vagnetti (2002) also constructed a 
similar spectral model to examine the color variability of 
AGNs in the sample of Giveon et al. (1999). 

We carried out multiwavelength monitoring obser- 
vations of many type 1 AGNs using the MAGNUM 
telescope with high photometric accuracy, and high- 
resolution imaging data obtained by HST are now avail- 
able for some of the MAGNUM targets. Therefore, us- 
ing the available high-resolution HST imaging data and 
the moderate-resolution ground-based images obtained 
by the MAGNUM telescope, we measure the flux contri- 
butions of the host galaxy for these AGNs with improved 
accuracy. Using previous studies, we also measure the 
flux contribution of narrow emission lines, which con- 
taminate the photometry and affect the color variability 



as well. Then, subtracting the combined contribution 
of the non- variable host galaxy plus narrow lines com- 
ponent from the optical flux of each target AGN, we in- 
tend in this study to settle a long-standing problem as to 
whether the spectral shape of AGN continuum emission 
in the optical region would change or remain constant 
with changing the brightness. 

In Section 2 we describe the characteristics of the tar- 
get AGNs, observations, reduction, and photometry. In 
Sections 3 and 4 we estimate the flux contributions of 
host galaxy and narrow lines, respectively. In Section 5 
we obtain the optical flux of the target AGNs and dis- 
cuss their color variability based on the linear correlation 
between the fluxes in any two different bands and the lo- 
cation of the non- variable component of host galaxy plus 
narrow lines in the flux to flux diagram. In Section 6 
we estimate the flux contribution of a dust torus to the 
observed optical flux. In Section 7 we discuss our inter- 
pretation of optical color variability of AGNs. In Section 
8 we summarize our results. 

2. MULTICOLOR LIGHT CURVE OF AGNS 

We describe the procedures of observation, data reduc- 
tion, and photometry only briefly. More details are given 
in Suganuma et al. (2006). 

2.1. Observations 

Our targets were selected from a sample of AGNs in- 
cluded in the MAGNUM program (Yoshii 2002; Yoshii, 
Kobayashi, & Minezaki 2003), on the condition that 
high- resolution HST/ ACS images were available and 
their host galaxy contamination could be accurately esti- 
mated. We used archival HST/ACS image data obtained 
by Bradley Peterson (Proposal ID: 9851, 10516). As a 
result, our targets consist of nine Seyfert 1 galaxies, one 
radio-quiet QSO, and one radio-loud QSO. All of them 
are nearby (z < 0.06) and are relatively low in luminosity 
{Mb > —22). Their characteristics are listed in Table 1. 

Monitoring observations were made using the multi- 
color imaging photometer (MIP) mounted on the MAG- 
NUM telescope (Kobayashi et al. 1998a, 1998b). The 
MIP has a field of view of 1'.5 x 1'.5 and obtains images 
simultaneously in the optical (UBVRI) and the near- 
infrared {JHK) bands by splitting the incident beam 
into two different detectors, an SITe CCD (1024 x 1024 
pixels, 0.277 arcsec pixel"^) and an SBRC InSb array 
(256 X 256 pixels, 0.346 arcsec pixel^^). 

The observational parameters of the target AGNs are 
listed in Table 2. Our targets were observed mostly in 
the U BVIJHK bands for seven years, but some of them 
were observed in fewer bands and for shorter periods. 
The typical sampling interval in the V band was 3 to 
30 days (median value). Observations were made most 
frequently for the (V, K) band pair: the second most fre- 
quently observed pair is {B,H), and the third is {U, J) 
or (/, J). This ordering was selected because of our aim 
at studying dust reverberation of AGNs as the first pri- 
ority (Minezaki et al. 2004; Suganuma et al. 2006). The 
FWHMs of the point-spread function (PSF) during the 
observations were typically 1".0-1".5 in the optical and 
0".8-l".0 in the near-infrared. Monitoring observations 
were performed with our own highly automated observ- 
ing system (Kobayashi et al. 2003). 
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2.2. Reduction and Photometry 

The images were reduced using IRAF. We followed the 

standard procedures for image reduction, with small cor- 
rections for nonlinear detector response. The nucleus 
fluxes of our targets were measured relative to the nearby 
reference stars, and then the fluxes of the reference stars 
were calibrated. Aperture photometry within (f) = 8". 3 
with the sky reference of a </> = 11". 1-13". 9 annulus was 
applied to all images, and then the fluxes of the nu- 
clei and reference stars were compared. We note that 
the sky annulus would include a significant portion of 
host-galaxy starlight for nearby target AGNs because 
of the apparently large extents of their host galaxies. 
The fluxes of the reference stars were calibrated with 
respect to photometric standard stars taken from Lan- 
dolt (1992) and Hunt et al. (1998). The uncertainty in 
the measurements was determined as follows: The rms 
variation of all measurements or the average photometric 
error for each image estimated from shot noise of pho- 
tons, sky background noise, and readout noise, whichever 
was larger, was divided by sqiiarc root of the number of 
images taken. Corrections for Galactic extinction were 
applied using the values given in the NASA/IPAC Extra- 
galactic Database (NED; Schlcgel, Finkbeiner, & Davis 
1998). Light curves based on aperture fluxes of all targcit 
AGNs are presented in Figures 1 to 11, where the offset 
fluxes of the host galaxy and narrow lines have not yet 
been subtracted. 

3. CONTRIBUTION OF HOST GALAXY 

In order to estimate the flux contribution of the host 
galaxy, we applied the surface brightness fitting using the 
two-dimensional image decomposition program GALFIT 
(Peng et al. 2002), which could fit the target image to 
analytic functions for the galaxy component of bulge and 
disk, plus a PSF for the nucleus and the sky contribution. 
Next, we applied aperture photometry with the same 
aperture size and the same sky annulus to a nucleus-free 
target image made by subtracting the best-fit PSF com- 
poncint from the target image. Since the surface bright- 
ness fitting technique would decompose the image into 
the galaxy component (in particular the bulge) and the 
nucleus component better when applied to higher spa- 
tial resolution images, we used HST/ACS high- resolution 
channel (HRC) data from the systematic AGN imaging 
survey by Bradley Peterson (Proposal ID: 9851, 10516). 
They were obtained in the F550M band. However, the 
flux contribution of the host galaxy in multiband is 
needed to be estimated to examine the color variation of 
the AGN optical continuum during flux variation. Thus, 
we applied the surface brightness fltting to the multiband 
MAGNUM images, fixing the host-galaxy shape parame- 
ters to those that were derived from the HST image, and 
applied aperture photometry to the nucleus-free MAG- 
NUM images. In this way, we expect to improve the ac- 
curacy of host-galaxy estimation in ground-based images 
of lower spatial resolution. 

3.1. Decomposition of HST/ACS/HRC Images into 
AGN and Host Galaxy 

Although Bentz et al. (2006, 2009) applied the sur- 
face brightness fitting to the HST imaging data to decom- 
pose them into two components of AGN and host galaxy, 



we here perform such fitting independently and compare 
both results in Section 3.3. The archival HST data con- 
sist of three images with different exposure times (120 s, 
300 s, 600 s) for each object. Since some of them were sat- 
urated, we used the unsaturated image with maximum 
exposure time for each object.'' The image identifica- 
tions for target AGNs arc fisted in Table 3. The PSF 
image was generated using the Tiny Tim package (Krist 
1993), which could model the HST optics plus the camera 
and filter system specifications. The input parameters of 
Tiny Tim, the location on HRC and the spectral shape 
of the incident beam, were set to the nucleus location 
and the power-law spectrum of fi, oc i/° (Wilhite et al. 
2005; Tomita 2005; Tomita et al. 2006; Suganuma et 
al. 2006), respectively. 

The HST images were corrected for the geometric dis- 
tortion using the IRAF task of PyDrizzlc. First, the 
surface brightness fitting was applied to the distortion- 
uncorrected images with the distortion-uncorrected PSF 
image generated by Tiny Tim and the analytic functions 
for bulge and disk of galaxy model. This procedure sig- 
nificantly decreased the fitting rcsidiial around the nu- 
cleus, compared to the fitting applied to the distortion- 
corrected HST image with the distortion-corrected PSF 
image. Then, the correction for the distortion was ap- 
plied to the best-fit model-galaxy image. Finally, surface 
brightness fitting was applied to the distortion-corrected 
model-galaxy image to derive the best-fit parameters of 
galaxy model (the fitting residuals here were found to be 
extremely small). 

The galaxy component was decomposed into bulge and 
disk fitted to a dc Vaucouleurs' (1948) R^''^ profile and 
an exponential profile, respectively, while only the de 
Vaucouleurs' profile was used for Mrk 110, Mrk 590 and 
PG0844-I-349, because the decomposition into bulge and 
disk did not give a convergent result with any initial val- 
ues of the parameters. We performed the galaxy de- 
composition by setting almost all parameters to be fitted 
except for the disk parameters; Since the disk param- 
eters often converged at unrealistic values or diverged 
without any restriction, we fixed the disk center at the 
nucleus location and constrained the disk scale length 
and axis ratio within the respective ranges that were de- 
fined from the major and minor isophotal diameters of 
the target galaxy taken from NED. However, we adopted 
the best-fit disk scale length to be much smaller than 
those expected from the isophotal diameter for NGC 
3516, NGC 4593, and 3C120, because small-scale struc- 
ture in the central region with complicated profile seemed 
to be practically fitted by an exponential profile for NGC 
3516 and NGC 4593 (Ferruit et al. 1998; Kormendy et 
al. 2006), and because the isophotal diameter of a bulge- 
dominated galaxy would not be a good indicator for the 
disk scale length for 3C120. Initial values of the parame- 
ters were set at those roughly estimated with IRAF, and 
then slightly different values were selected to examine 

We also applied the surface brightness fitting to the HST imag- 
ing data with three different exposure times by stacliing them fol- 
lowing Bentz et al. (2006), and found little difference in the model 
parameters determined for the decomposed components. In ad- 
dition, the resulting host-galaxy fluxes using our decomposition 
parameters and those using the parameters derived by Bentz et 
al. (2009) are found to be consistent with each other, as will be 
described in Section 3.3. 
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how final values of the parameters varied with different 
initial values. The sky level of the parameters was fitted 
with its initial value set to be zero. When the fitted val- 
ues were unrealistic, the best-fit values were searched for 
by changing the sky level manually, until the resultant 
sky fluxes were found to be comparable with those calcu- 
lated from the HRC user's manual. Finally, we adopted 
the values of parameters that minimized the reduced x^- 
Figures 12 to 14 show the best-fit and residual images 
of the HST data. The best-fit parameters of the galaxy 
components are listed in Table 4. 



3.2. Decomposition of MAGNUM Images into AGN 

and Host Galaxy 

To accurately estimate the host galaxy, we selected a 
moderate number of observations for each target and fil- 
ter from MAGNUM monitoring data taken under good 
and stable seeing conditions and in the AGN-faint phase. 
The PSF image was made from the images of nearby ref- 
erence stars for each target AGN. These stellar images 
were not taken simultaneously with the AGN images, 
so stable seeing was required for accurate fitting. In this 
way, more than five nights of images were selected in each 
band, although for some AGNs fewer than five nights 
were selected in each band because of the lack of appro- 
priate data. Then, we combined the AGN images and 
stellar images from each night, respectively, to increase 
the signal to noise ratio and clipped the field region that 
is the same as the HST field for each AGN. 

These combined and clipped images of target AGNs for 
each night were fitted to the PSF made as described for 
the nucleus component and to analytic functions for the 
galaxy component of bulge and disk, with the parameters 
of host-galaxy shape (scale length and axis ratio) fixed at 
the best-fit values of the surface brightness fitting using 
the HST data. The sky level was fixed at the mean sky 
value of reference star images. When was obviously 
large, we did not fix the disk scale length. This occurred 
mainly in the B and / images, probably because of pos- 
sible radial color gradient in a galaxy, arising either from 
stellar population gradient or from metallicity gradient 
(e.g., de Jong 1996; Taylor et al. 2005). However, this 
had eventually little eflFect in estimating the host galaxy 
flux (differences were less than a few percent). Finally, 
we applied the aperture photometry with the same aper- 
ture size to a nucleus- free image made by subtracting the 
best-fit PSF component from the original image. With 
a moderate number of independent observations avail- 
able, the host-galaxy flux was estimated as the average 
of those fluxes of all observations, and its error was esti- 
mated as their standard deviation. Figures 15 to 17 show 
the best-fit and residual images of the MAGNUM data. 
The estimated values of the host-galaxy flux for target 
AGNs within the monitoring aperture are listed in Table 
5. 



3.3. Consistency Check of Host Galaxy Estimates 

In order to improve the estimation of the host-galaxy 

flux in the ground-based image, we decomposed the 
MAGNUM image into the components of AGN and host 
galaxy using the shape parameters of host galaxy derived 



from the HST image. Here, we compare our results of 
the host-galaxy flux from the MAGNUM images in the V 
band with those estimated from the HST images in the 
F550M band, because the effective wavelengths of these 
bands are very close to each other. We note that they are 
based on the aperture photometry with the same aper- 
ture size to the nucleus-free image made by subtracting 
the best-fit PSF component from the HST image. The 
procedure of image decomposition is described in Section 
3.1. The error of the host-galaxy flux was estimated as 
follows: An artificial image was generated with the best- 
fit parameters by adding artificial noise. Then, the PSF 
flux of the artificial image was increased and decreased 
with the total fiux of PSF and bulge kept constant, and 
the reduced was calculated each time. Finally, the er- 
ror was estimated as the range of the host-galaxy flux of 
the artificial image such that the reduced x^ calculated 
above is the same as that of the observed image fitted 
with the best-fit parameter values. 

Figure 18 shows the host-galaxy fluxes of the 11 target 
AGNs in the V band for MAGNUM and in the F550M 
band for HST. We have found that both results are con- 
sistent with each other within a 3(t level of errors. We 
have also found that the systematic difference between 
them is as small as 1.4 % and the standard deviation is 
18 %, which further decreases to 11 % when the eight 
targets with brighter host galaxies are selected. This 
good agreement suggests that our estimation of the host- 
galaxy flux in the MAGNUM image provides an accuracy 
comparable to that in the HST image. 

Bentz et al. (2009) performed the surface brightness 
fitting of host galaxy component using the Sersic (1968) 
profile. Since his method allows a more general distribu- 
tion compared with our use of de Vaucoulcurs' profile and 
an exponential profile, it is instructive to decompose the 
MAGNUM image in the same way as described in Sec- 
tion 3.2, but with the parameter values of host-galaxy 
shape in Table 4 in Bentz et al. (2009). Figure 19 shows 
the host-galaxy fluxes of 11 target AGNs with the pa- 
rameter values of host-galaxy shape adopted by Bentz et 
al. (2009), in comparison of those with our own values of 
such parameters. We have found that the systematic dif- 
ference is only 0.3% and the standard deviation is 4.3% 
for the eight targets with brighter host galaxies, while 
for other three targets with fainter host galaxies, both 
results are consistent with each other within a 2a level 
of errors. This good agreement between the host-galaxy 
fluxes estimated with two sets of host-galaxy parameters 
suggests that the uncertainty of AGN and galaxy de- 
composition in the HST images is not significant in our 
analysis below. 

3.4. Host Galaxy Color 

The estimated colors of host galaxies within the pho- 
tometric aperture are listed in Table 6. For compari- 
son, typical colors of galaxies (Fukugita et al. 1995) and 
bulge (Kinney et al. 1996) are also presented in Ta- 
ble 6. For nearby targets (NGC 4051, NGC 3227, NGC 
3516, NGC 4593, NGC 5548, Mrk 590, Mrk 817), the 

* Since wc intended to estimate tlie host-galaxy flux within the 
aperture around the center and not to study its whole structure, 
we did not fit the high-resolution HST image and the wide-area 
MAGNUM image together. 
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dominant host-galaxy component within the aperture is 
identified with a bulge, and the B — V and V — I col- 
ors are comparable to those colors of bulge or slightly 
bluer. This result is interpreted as meaning that the 
host flux within the aperture is mainly from the bulge 
or slightly contaminated by the disk, which is generally 
bluer than the bulge. For NGC 4151, the B - V color 
is interpreted as above, but the V — I color is slightly 
bluer than the bulge. However, the small contamina- 
tion of the extended [O III] emission into the host-galaxy 
fluxes in the B and V bands could explain the colors of 
NGC 4151 's host galaxy. Observation by Yoshida et al. 
(1993) show that the [O III]A4959, 5007 emission lines of 
NGC 4151 come not only from the nucleus but also from 
the extended region (~10" x 10" around the nucleus), and 
the flux of extended [O III] emission is about 10-20% of 
nucleus [O III] emission. The [O III] emission that con- 
tributed to the y-band flux would make the V — I color 
slightly bluer, while it happens to contribute to both the 
i?-band and V-band fluxes similarly to the color of the 
host galaxy. Consequently, the B — V color stays almost 
unchanged. 

For relatively distant targets (3C 120, Mrk 110, PG 
0844-1-349), the flux of the host galaxy within the aper- 
ture contains nearly all the flux from the host galaxy. For 
Mrk 110, the dominant host-galaxy component within 
the aperture is identified as a late-type spiral disk, and 
the B — V and V — I colors are comparable to those of a 
typical spiral galaxy, but for 3C 120 and PG 0844+349, 
although these dominant host galaxies are identified as 
early-type galaxies, their color tends to be bluer than 
that of a typical early-type galaxy. However, many au- 
thors have reported that the host-galaxy colors of QSOs 
tend to be bluer than that of a typical early-type galaxy 
(Kauffmann et al. 2003; Sanchez et al. 2004; Schweitzer 
et al. 2006; Jahnke et al. 2004, 2007). Therefore, the 
host-galaxy colors we estimated are consistent with the 
trend reported by previous studies. 

4. CONTRIBUTION OF NARROW LINES 

Observations of nearby AGNs show that the spatial 
extent of the narrow-line region is typically around 100 
pc or larger (Schmitt et al. 2003), so the narrow-line flux 
does not respond to variations in incident flux from the 
nucleus and is generally regarded as non-variable over 
100 yr. Therefore, it is possible to estimate the narrow- 
line flux from past spectral data. We considered only 
strong narrow fines ([O III] A4959, 5007, H/3, H7), because 
these lines occupy a large fraction of the optical narrow- 
line flux (about 90% or more, except for the small blue 
bump) and the flux of other lines is hardly estimated. 

Compilation of past [O III]A5007 observations by 
Whittle et al. (1992) indicates that the narrow-line 
fluxes of very nearby AGNs within different apertures 
are not always the same. The scatter of narrow-line 
fluxes in large apertures of 7"-15" is smaller than that 
in small aperture of 2"-4". They interpreted this differ- 
ence as the seeing effect and angular size effect (a large 
aperture covers most of extended narrow line regions). 
For [O III]A5007, we adopted the mean fluxes of large- 
aperture observations presented in Whittle et al. (1992), 
because the aperture size is comparable to our photo- 
metric aperture {(p — 8". 3), and the scatter of narrow- 
line fluxes is small in many large-aperture observations. 



For Mrk 110, Mrk 590, and Mrk 817, we adopted the 
mean fluxes of small-aperture observations in Whittle et 
al. (1992) because large-aperture data were not avail- 
able. We note that the quality of these data is poorer. 
[O III]A5007 data for 3C 120 and PG 0844-h349 were 
not presented in Whittle et al. (1992), so we adopted 
data from Tadhunter et al. (1993) obtained with a 5" 
slit and from Miller et al. (1992) obtained with a 6" slit, 
respectively. 

For [O III]A4959, we adopted fluxes computed with a 
theoretical line ratio of [O III]A5007/[O III]A4959 = 3.01 
(Storey & Zcippcn 2000). For H/?, we adopted fluxes for 
some targets from Kaspi et al. (2005), and for others 
those computed with a roughly estimated line ratio of 
[O III]A5007/H/3 ~ 10 (Veilleux & Osterbrock 1987). For 
H'y, we adopted fluxes computed with a theoretical line 
ratio of H7/H/3 ^ 0.45 (Case B: low-density limit and 
T = 5000 K, Osterbrock & Ferland 2005). All fiuxes 
were corrected for Galactic extinction using the linearly- 
interpolated NED value. The adopted values of narrow- 
line flux are listed in Table 7. 

Finally, we estimated the flux contribution of narrow 
lines in the wavelength coverage of each filter by con- 
volving the narrow-line spectrum with the MAGNUM 
filter transmission curve. We assumed a 5-function for 
narrow-line profiles for simplicity, because the width of 
the narrow lines is so small that the filter transmission 
can be taken as constant at their wavelengths. The esti- 
mated values of fiux contribution of the narrow lines in 
the B and V bands are listed in Table 8. We note that 
the narrow-line component has a minor contribution to 
the total flux in aperture, and the host galaxy component 
always makes a major contribution. 

5. EXAMINATION OF COLOR VARIABILITY 

In order to examine the color variability of optical con- 
tinuum emission of type 1 AGNs with their flux varia- 
tion, we applied flux to flux plot analysis, originally pro- 
posed by Cholonicwski (1981). Using the UBV monitor- 
ing data of Seyfert galaxies, Cholonicwski (1981) plotted 
the data of each target obtained at different epochs on 
a diagram of the total flux in one band versus the total 
flux in another band, and first noticed a linear correla- 
tion between the fluxes of an AGN in two different optical 
bands. This has subsequently been confirmed by many 
authors (Winkler et al. 1992; Winkler 1997; Suganuma 
et al. 2006; Tomita et al. 2006), and it is assumed that 
the vector describing the total observed flux in the flux 
to flux diagram is the sum of a constant flux vector of 
starlight of the host galaxy and a variable flux vector 
of the AGN continuum that is constant in direction but 
not in magnitude. Assuming a variable flux vector of 
the AGN continuum that does not change its direction 
is equivalent to assuming the constant spectral shape of 
the AGN continuum. 

We first examine the linear correlation between the 
fluxes in two different bands by plotting our long-term 
monitoring data with high photometric accuracy in a flux 
to flux diagram. Then, by siunniing the host-galaxy flux 
and narrow-line contribution, we locate the non-variable 
component in the flux to flux diagram. In this way, we 
are able to verify the simple two- vector model, by exam- 
ining whether the non-variable component is located on 
the fainter extension of the linear correlation line. 
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5.1. Flux to Flux Plot and Linear Fit 

In this study, we plot only {B, V) and {V, I) pair data 
in the flux to flux diagram, because our monitoring data 
are largest in the V band, and the number of other paired 
data obtained on the same nights is relatively small. 
J7-band data were not used because it was difflcult to 
estimate the significant contribution of the small blue 
bump, blended Fe II lines and Balmer continuum in this 
band. Figures 20 and 21 show the B-band flux to V- 
band flux diagrams for all 11 targets. Figures 22 and 23 
show the F-band flux to /-band flux diagrams for 8 tar- 
gets, because the other 3 targets, NGC 3516, NGC 4593, 
and PG0844-I-349, were not monitored in the / band. 
Straight-line fitting to the pair data in the fiux to fiux 
diagram was done with the fitting code made by Tomita 
(2005) following the most generalized multivariate least 
square process as given by Jefferys (1980, 1981). The pa- 
rameter values of the best-fit linear regression lines, the 
number of the data points, and the reduced are listed 
in Table 9. The residuals of the data from the best-fit lin- 
ear regression line are plotted in the lower panels of the 
flux to flux diagrams in Figures 20 to 23, and the ratio of 
Ry = cr^^VAF is fisted in Table 9, where a'^^ is the stan- 
dard deviation of the residuals and AV = /y — /y"" is 
the amplitude of fiux variation. 

As shown in Figures 20 and 21, the {B,V) pair data 
are certainly distributed in a highly linear manner for all 
11 targets, and no curvature can be seen in the plots. As 
listed in Table 9, the residuals of the straight-line fitting 
of the {B, V) plots are so small that the reduced is 
near unity, except for two targets of NGC 4151 and Mrk 
110. Indeed, according to the x^-test for the straight- 
line fitting of the {B, V) plots, the linear relationship is 
not rejected for seven of the 11 targets at a 1% level of 
significance. On the other hand, for the rest of the tar- 
gets, NGC 4051, NGC 4151, NGC 5548, and Mrk 110, 
the x^-test rejects the linear relationship at the 1% level 
of significance, but the ratio of fitting residual to flux 
variation Ry is comparable to or even less than that for 
the other seven targets. The x^-test and the i?y-ratio 
suggest that higher photometric accuracy and/or larger 
number of data could reveal small deviation from the lin- 
ear correlation for these objects. We here note that small 
contamination by broad emission- lines, which show time- 
delayed response to optical continuum emission, is liable 
to cause a small scatter around the linear correlation for 
the (B, V) plots. 

As shown in Figures 22 and 23, the {V,I) pair data 
are also certainly distributed in a highly linear manner 
for all eight targets, and no curvature can be seen in 
the plots. However, the reduced x^-value of the straight- 
line fitting and the i?v-ratio of fitting residual to flux 
variation for the {V, I) plots tend to be larger than those 
for the {B, V) plots. Indeed, according to the x^-test, 
the linear relationship for the (V, /) plots is rejected for 
five of eight targets at the 1% level of significance. The 
origin of the relatively larger scatter in the {V, I) plots 
will be discussed in Section 6. 

5.2. Location of Host- Galaxy plus Narrow-Line 
Component in the Flux to Flux Diagram 

The location of the host-galaxy plus narrow-line com- 
ponent (HOST-I-NL component) in the flux to flux di- 



agram is presented in Figures 20 to 23 for the target 
AGNs. We define th e error of the HOST+NL flux as 
a = V^host+^'nl' where cthost is the error of the 
host-galaxy flux and ctnl is the error of the narrow-line 
flux, respectively. We assume 10% for ctnl, because the 
error of the narrow-line flux is not presented in Whittle 
(1992) but the standard deviation of several measure- 
ments of the narrow-line flux for the individual objects 
in his sample ranges from several % to a few of 10 %. 
We note that ctnl usually has a minor contribution to 
the total error a. It is evident that the HOST-I-NL com- 
ponent is located on the fainter extension of the best-fit 
linear regression line within 2a error in almost all cases, 
and even in only one exception of Mrk 817 (V,/), the 
deviation of the HOST-I-NL component from the best- 
fit linear regression line is still about 3a level. We note 
that the angular size of bulge component of Mrk 817 is 
so small {Re = 0".5) that decomposing the bulge and 
nucleus is relatively diflacult. Therefore, the linear cor- 
relation between the fluxes in different bands and the 
location of the HOST-j-NL component on the flux to flux 
diagram strongly indicate that the AGN optical contin- 
uum does not systematically change the spectral shape 
while changing the brightness. In other words, it does 
not become bluer even in the brighter phase. 

By subtracting the HOST-j-NL component, we arrive 
at a conclusion that during seven years of monitoring 
observations, the AGN optical contimium fiux varied by 
0.3 3.0 mag. In fact, the optical spectral shape of the 
AGN continuum emission of NGC 4151 and NGC 5548, 
which showed the largest-amplitude variations, remained 
nearly constant during flux variations of more than 10 
times. 

6. EFFECTS OF DUST TORUS EMISSION 

While a straight line provides a good fit to almost all 
{B, V) and (V, I) plots, we see noticeable scatter around 
the line for many of the (V, I) plots. This was also 
pointed out by Winkler (1997), who suspected the exis- 
tence of some additional variable component that could 
contaminate the 7-band flux and cause scatter in the 
{V, I) plots. 

Thermal emission from dust torus is thought to be the 
most promising candidate for an additional variable com- 
ponent that could contaminate the optical flux at longer 
wavelengths. This emission in the / band could then dis- 
turb the linear relationship in the (V, /) plots because it 
shows a time-delayed response to emission in the same 
band from the accretion disk (Minczaki et al. 2004; Sug- 
anuma et al. 2006). We here estimate the flux contribu- 
tion of dust torus in the / band, using simultaneous VIJ 
light curves obtained by the MAGNUM telescope. 

First, we assume that the flux in each band is decom- 
posed into accretion disk, HOST+NL, and dust torus, 
while the emission from dust torus in the V band is 
ignored. That is because the sublimation temperature 
of dust is about 1800 K and docs not exceed 2000 K 
(Salpetcr 1977; Fodders & Fegley 1999), and blackbody 
radiation at 2000 K in the V band is less than 5% of that 
in the / band. Then, the fluxes in the VLJ bands are 
expressed as 

fv{t) = 0''{t) + fr , (1) 

fi{t)=ff''{t)+fr+fP^{t) , (2) 
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where the superscript AD stands for accretion disk, HN 
for HOST+NL, and DT for dust torus. Second, we as- 
sume that the V to J fliix ratio in the power-law spectral 
emission from the accretion disk is kept nearly constant 
irrespective of its flux variation, namely 



fntyfnt) 



(4) 



where vy and i>j are the effective frequencies in the 
MAGNUM V and J filters, respectively, and ayj is the 
power-law spectral index of constant value. Finally, al- 
though the constant color of / — J and the synchronized 
variation in the I and J bands for dust torus are still un- 
der discussion, we assume that the I to J flux ratio in the 
thermal emission from dust torus is also kept constant, 
namely 



f?^{t)/fy-{t) 



DT/ 
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where the power index of j3 has a constant value of 1.6 
from the absorption efficiency Qabs for graphite parti- 
cles of radius 0.05/xm in the relevant wavelength region 
(Draine & Lee 1984), and B^{Td) is the specific inten- 
sity of a black body of constant dust temperature Ta 
at frequency v. We note that a single temperature is 
assigned to dust torus, which is realistically of a multi- 
temperature component. 

With the above assumptions, the /-band flux from dust 
torus is written as 

fF^it) = ( fj{t) - ifvit) - f^''){yj/yvr-^ - fD 

x(vilvjf{B,ATd)IB,ATd)) , (6) 

where Jv{t) and fj{t) are the observed light curves, fy^ 
was estimated in Section 5.2, and was estimated as 
described in Section 3 using MAGNUM J-band images 
(the narrow-line contribution in the J band was ignored) . 
Consequently, two parameters, and ayj, remain to be 
determined. In other words, given these parameters, we 
can subtract the contribution of dust torus fp^{t) from 
the observed /-band flux and fit a straight-line to the 
(V, /) plots corrected for dust torus emission. 

We searched for the best set of these parameters that 
could minimize the x^-value for thc^ straight-line fitting 
to such dust-free (V, /) plots by repeating the calcu- 
lation from Td = 1000 K to 3000 K at intervals of 50 K 
and from ayj = —3 to -|-3 at intervals of 0.1. For exam- 
ple, Figure 24 shows the contour map of reduced for 
the parameter search for NGC 4151. Figures 25 and 26 
present the y-band flux to /-band flux diagram and the 
best-fit linear regression line with and without the dust 
correction for seven targets. Mrk 590 was removed from 
the discussion because its AGN component was so much 
fainter than the host galaxy in the J band that accurate 
image decomposition was difficult. The residuals of the 
data from the best-fit linear regression line are also plot- 
ted in the middle and lower panels of Figures 25 and 26. 
The best-fit parameters of and ayj with the J-band 
host-galaxy flux f^^ are tabulated in Table 10, together 
with the best-fit ayi and the reduced of the straight- 
line fitting to the {V, I) plots after subtracting the dust 



component in the /-band. For reference, the best-fit ayi 
and the reduced without the dust correction are also 
tabulated in Table 10. 

As a result of the correction for contamination by dust 
torus emission in the / band, we found that the re- 
duced x^-values for all target AGNs dramatically de- 
creased to x^ ~ 1 foi' six of the seven targets except 
for NGC 4151. Indeed, the linear relationship for the 
{V, I) plots after the dust correction turned out to be 
not rejected at the 1% level of significance, and even for 
NGC 4151, the reduced x^-value decreased remarkably 
with the dust correction. We also found that the location 
of the HOST-I-NL component in the {V, I) plot became 
closer to the fainter extension of the best-fit regression 
line after the dust correction in almost all cases. In ad- 
dition, the best-fit value of Td (1700 ~ 2000 K) became 
consistent with the \ipper limit of the sublimation tem- 
perature of dust grains. On the other hand, the best 
value of avj (—0.7 ~ 0.0) was generally smaller than 
the value reported by Tomita (2005). Most values of 
ayj showed good agreement with the value of ayi, but 
not for NGC 5548 and Mrk 817. However, ayj woifld 
not be constrained well because the reduced x^-value was 
insensitive to any change of ayj- The reduced x^-value 
for Mrk817 was so small that the linear relation seemed 
to provide a satisfactory fit without the dust correction 
®, however, it also much decreased by the dust correc- 
tion, and the best-fit Td showed a reasonable value. This 
suggests that the /-band fiux is contaminated by a dust 
component for Mrk817 as well. 

Prom these exercises, we conclude that the contribu- 
tion of dust torus to the optical fiux at longer wavelengths 
cannot be negligible in disciissing the optical color vari- 
ability, and most of the scatter around the best-fit re- 
gression line in the {V, /) plot can be explained by con- 
tamination of dust torus emission to the observed /-band 
fiux. Along with the results of Section 5, this result con- 
firms that the optical spectral shape (0.44 ~ 0.79 i^m) 
of AGNs remains nearly constant during their flux vari- 
ation. 

7. DISCUSSION 

7.1. Spectral Variability of UV- Optical Continuum 
Emission in Previous Studies 

The results of preceding sections showed that AGN op- 
tical continuum does not systematically change in spec- 
tral shape during flux variation. In this section, we re- 
view the previous studies of AGN spectral variability in 
the UV-optical region and discuss them in comparison 
with our results. 

7.1.1. Optical Color Variability by Webb & Malkan 
(2000) and Giveon et al. (1999) 

We examine two leading studies that claimed spectral 
hardening of intrinsic emission from AGNs, opposite to 
our conclusion. 

Webb & Malkan (2000) monitored 23 AGNs in the 
blue, yellow, and red bands for three months. They con- 
structed a spectral variation model of AGN plus host 
galaxy with the AGN component contributing two thirds 

^ We note a small number of the (V, /) points available for Mrk 
817. 
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of the total flux in the yellow band as an extreme case for 
a large contribution of host galaxy. Then, allowing 0.5 
mag variation for the AGN component, they estimated 
the amplitude of flux variation in the blue, yellow and 
red bands, and found, as in their Figure 8, that the blue 
to red amplitude ratio exceeded that expected from the 
spectral variation model for two of six targets. Accord- 
ingly, they concluded that the intrinsic emission from 
AGNs became bluer when they were brighter. 

Webb & Malkan (2000) estimated the host-galaxy con- 
tribution for more than two thirds of their objects. Since 
four of them are included in our sample, it is instructive 
to compare their estimation of host-galaxy flux with our 
estimation with improved accuracy as described in Sec- 
tion 3.3. Their host-galaxy flux in the red band is about 
13 mJy (NGC 3227), 8.3 mJy (NGC 4151), 3.0 mJy 
(NGC 5548), and 1.4 mJy (Mrk 817), with an absolute 
calibration error of 20 — 30%. Their aperture diameter 
of 11" is larger than ours of 8". 3, but their estimation 
of host-galaxy flux is much smaller than ours in the / 
band (the effective wavelength of their red band is close 
to that of the / band). Moreover, their assumption of 
the host-galaxy component contributing one third of the 
total flux in the yellow band is not at all an extreme case. 
Table 11 shows the contribution of the HOST-I-NL com- 
ponent in the average flux of our light curve data. We 
sec that the contribution in the V band, for which the 
effective wavelength is slightly shorter than that of the 
yellow band, exceeds one third for eight out of 11 targets, 
and even two thirds for 5 targets. Therefore, the targets, 
which show a large blue to red variation amplitude ratio 
in Webb & Malkan (2000), do not necessarily indicate 
the spectral hardening of intrinsic emission from AGNs. 

Giveon et al. (1999) monitored 42 PG quasars in the 
BR bands for seven years, which were generally more 
luminous than our targets. Their photometry was based 
on PSF fitting, and the host-galaxy contribution was not 
subtracted. They found that at least half of their targets 
showed bluer B — R color when they became brighter. It 
is therefore worthwhile to examine the color variability 
of their targets by means of linear correlation analysis 
described in Section 5. 

First, for PG 0844+349 that is the only target in com- 
mon with ours, a tight linear relationship for the (B, V) 
plot is presented and the HOST+NL component is found 
to be located on the fainter extension of the best-fit lin- 
ear regression line. Second, for six targets that were 
more luminous than PG 0844+349 by about 1.5 mag 
in the B band and that showed significant flux varia- 
tion of AB > 0.6 mag (see their Table 3), their BR 
monitoring data are plotted in a flux to flux diagram 
as shown in Figure 27, together with the best-flt regres- 
sion line. We sec that the linear relationship is a good 
approximation for the six targets, in disagreement with 
the conclusion by Giveon et al. (1999) that they all 
became bluer when they became brighter. The scatter 
around the best-fit regression line in Figure 27 is some- 
what larger than that for our targets shown in Figures 
20 to 23. Such scatter would be caused by the broad 
emission-line components of small blue bump and Balmer 
lines that contaminate the flux in the B and R bands, 
and also by the uncertainty associated with PSF-fitting 
photometry. The good linear relationship for the {B, R) 
plots demonstrates that the optical spectral shape of the 



variable component remains nearly constant even for the 
luminous QSOs. Although it is not examined whether 
or not the HOST+NL component of them locates on the 
fainter extension of the linear regression line, by combin- 
ing our results for Seyfert galaxies and PG 0844+349, 
it is strongly indicated that optical continuum of AGNs 
from Seyfert galaxies to QSOs docs not systematically 
change in the spectral shape during fiux variation. 

7.1.2. UV Color Variability 

Seyfert 1 galaxies and QSOs produce strong contin- 
uum emission in the UV region as well as in the optical 
region. Similarly to the optical color variability, two op- 
posite claims on the UV color variability have prevailed 
without reaching a conclusion as to whether or not the 
UV contimuim of AGNs changes in spectral shape during 
flux variation. 

Paltani & Walter (1996) observed the UV spectra of 
15 nearby AGNs at different epochs using the lUE satel- 
lite. By principal component analysis, they explained 
UV variation in almost all their targets as the sum of 
a variable component of constant power-law spectral 
shape, with evidence of reddening by dust extinction, 
and a non-variable component of small blue bump con- 
sisting of steep Balmer continuum and Fe II pseudo- 
continuum. Other studies arriving at similar results in- 
clude Rodriguez-Pascual et al. (1997) and Santos-Lleo 
et al. (1995), who observed Fairall 9 and NGC 4593, 
respectively, using the lUE satellite. From discussions 
on the linearity of UV monitoring data in the /(1803A) 
to /(1447A) diagram, they concluded the constant spec- 
tral shape in the UV region during their flux variation. 
Santos-Lleo et al. (1995) also estimated the non- variable 
component for NGC 4593 and found that it is located on 
the fainter extension of the best-flt regression line within 
the uncertainties. 

On the other hand, for NGC 5548 (Wamsteker et al. 
1990; Korista et al. 1995), NGC 4151 (Crenshaw et al. 
1996), and NGC 3783 (Reichert et al. 1994), the claim of 
spectral hardening in the UV region was argued from the 
changeover of the power-law spectral index through flux 
variation, or from the larger amplitude of variability at 
shorter wavelengths. Although most of these studies did 
not take into account the effect of non-variable compo- 
nent, Wamsteker et al. (1990) considered the host-galaxy 
component. Romano & Peterson (1998) plotted quasi- 
simultaneous UV and optical monitoring data for NGC 
5548 in the F(1350l) to i^(5100l) diagram. They found 
a positive correlation with a curvature to be statistically 
consistent with a linear correlation in the F {1350 A) - 
i^(5100A) plots, hence they argued for a change in spec- 
tral shape. 

Vanden Berk et al. (2004) studied the relationship 
between variability amplitude and rest-frame wavelength 
of the UV-optical region from SDSS two-epoch multicolor 

observations of about 25,000 quasars and found the larger 
amplitude at shorter wavelengths in the UV region of A < 
4000j4, although this trend seemed to be very small, and 
the amplitude was almost constant in the optical region 
of 4000A < A < 6000 A, as presented in their Figure 
13. In their study, the non-variable component was not 
subtracted, but they noted that if its contribution had 
been signiflcant, the variability amplitude should have 
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decreased drastically at A > 4000A, which would simply 
strengthen the correlation they found. Therefore, they 
concluded that quasars were systematically bluer when 
brighter. On the other hand, their results in the optical 
region of 4000A < A < 6OOOA, which overlaps with our 
B-band and V-band wavelength regions, show little or 
no change in amplitude, and in fact, they arc consistent 
with our conclusion of nearly constant spectral shape in 
the optical region. 

Although a problem is still controversial as to whether 
or not the UV spectral shape really changes, it is possi- 
ble to conceive a theoretical model that allows spectral 
hardening in the UV region and constant spectral shape 
in the optical region. Recently, Pereyra et al. (2006) 
presented such a model, assuming the standard accre- 
tion disk (Shakura & Syunyaev 1973) and variation of 
mass accretion rate onto the disk. The variation of accre- 
tion rate changes the maximum disk temperature with- 
out changing the radial temperature profile in the disk. 
Thus, the UV spectrum would become harder when it 
brightens, because it originates from the innermost disk 
where the maximum temperatiire increases with the mass 
accretion rate. On the other hand, the optical spectrum 
would remain nearly constant, because it originates from 
relatively outer region of the accretion disk where the 
temperature structure remains unchanged. Figure 1 of 
Pereyra et al. (2006) shows that the model predicts large 
spectral variation in the UV region while little spectral 
variation in the optical region. 

In addition, the color variability in the UV region 
might reflect different physical properties of individual 
AGNs. UV spectral hardening was reported for some 
AGNs (NGC 5548, NGC 4151, NGC 3783), while nearly 
constant spectral shape was reported for other AGNs 
(Fairall 9, NGC 4593). According to the model for flux 
variation above, the spectral shape of the UV contin- 
uum changes with the maximum disk temperature, which 
is dependent on the mass accretion rate and the cen- 
tral black hole mass. In order to settle this problem, 
highly accurate simultaneous UV and optical monitor- 
ing observations are necessary for many AGNs, together 
with high-resolution images for estimation of their non- 
variable component. 

7.2. Optical Color of Accretion Disk 

The slope of the best-fit regression line in the flux to 

flux diagram corresponds to the color of variable AGN 
component. Our conclusion in the preceding sections is 
that the optical spectral shape of AGNs remains nearly 
constant during flux variation, hence the optical color 
derived from the flux to flux diagram can be regarded as 
the color of the AGN optical continuum. Thus, the color 
of AGN optical continuum between the X and Y bands 
is derived as 

X - y = -2.5(logio(axy) - logio(,/xo//Yo)) , (7) 
where axY is the slope of the best-fit regression line in 
the X-h&nd flux to F-band flux diagram, and fxo and 
/yo si,re the fluxes at zero magnitude in the respective 
bands. Assuming the power-law spectrum of cx v" for 
the AGN optical contimium, we obtain the power-law 
spectral index a between the X and Y bands as 

log(axy) 



axY = 



(8) 



where i>x and i^y are the effective frequencies of the re- 
spective bands. 

Table 12 shows the color and spectral index obtained 
from the best-flt regression line of the flux to flux plots. 
They are compared with those obtained by Winkler 
(1997) and Vanden Berk et al. (2001), as shown in 
Figure 28. Furthermore, since the AGN optical contin- 
uum emission is considered to originate from the accre- 
tion disk (Shields 1978; Malkan & Sargent 1982), the 
spectral index expected from the standard accretion disk 
(Shakura & Syunyaev 1973) at long wavelengths with an 
infinitely extended disk is also shown in this figure for 
comparison. We found that the colors arc comparable 
to that of standard accretion disk or slightly redder, and 
are bluer than the QSO composite color by Vanden Berk 
et al. (2001). 

Winkler (1997) plotted the optical colors of the vari- 
able component of 92 Seyfert 1 galaxies in the two-color 
diagram, and found that most targets were distributed 
around the bluest target, while redder targets were dis- 
tributed in the direction of the reddening vector. He in- 
terpreted this result to mean that the bluest color was the 
intrinsic color of Seyfert galaxies, and the redder colors 
were caused by internal extinction from their circumnu- 
clear regions. The {B — V) to {V — I) diagram is shown 
in Figure 29, where the data of Winkler (1997) are plot- 
ted together with our data. As shown in Figure 29, our 
data present a similar distribution with those of Winkler 
(1997). Following his method, we estimated the internal 
extinction of redder targets and compared the extinction 
values with those estimated in previous studies. 

Paltani et al. (1996) estimated the internal extinction 
by fitting the power-law spectral emission to the UV con- 
tinuum for nearby Seyferts, and some of their targets 
were in common with ours. Crenshaw et al. (2001) esti- 
mated the internal extinction of NGC 3227 by comparing 
its UV spectrum with that of NGC 4151. We then com- 
pared our estimates of extinction with theirs and found 
them consistent with each other, as shown in Table 13. 
The consistency between the reddening of optical colors 
obtained from the flux to flux plots and the extinction es- 
timated from other methods supports the idea of Winkler 
(1997) that the optical color of the variable component, 
after corrected for internal extinction [E{B — V) ~ 0- 
0.2mag], is interpreted as the intrinsic color, which is 
similar to the color of standard accretion disk. 

8. SUMMARY 

We carried out accurate and frequent monitoring ob- 
servations of 11 nearby AGNs (9 Seyfert 1 galaxies and 
2 QSOs) in the B, V, and I bands for seven years with 
the MAGNUM telescope. We estimated the flux con- 
tribution of the non-variable HOST+NL component in 
the AGN optical continuum emission by applying surface 
brightness profile fitting to HST/ACS high-resolution im- 
ages and MAGNUM images for the host galaxy (HOST), 
and by analyzing the results of previous spectral obser- 
vations for narrow emission lines (NLs). We found a 
tight linear flux to flux relationship for the (B, V) and 
{V,I) plots observed on the same days for each tar- 
get AGN and also found that the HOST-I-NL compo- 
nent estimated above could be located on the fainter 
extension of the best-fit regression line in the flux to 
flux diagram within 2a uncertainties in almost all cases. 



10 



Sakata et al. 



These results strongly support the idea that AGN opti- 
cal continuum emission does not systematically change 
its spectral shape with changing flux. By subtracting 
the contribution of the HOST+NL component from the 
observed flux, we estimated the accurate variability am- 
plitude of AGN optical continuum emission and found 
that some AGNs (NGC 4151, NGC 5548) changed their 
optical brightness by more than 3 mag during the seven- 
year monitoring period. This means that some AGNs 
show optical flux variation by more than 10 times with 
their color kept nearly constant on a timescale of several 
years. Therefore, spectral hardening, which was argued 
otherwise in many studies, can be interpreted as an ap- 
parent trend due to contamination of the non-variable 
HOST+NL component in the optical region of 4400A- 
7900 A considered. The nearly constant optical color of 
the variable component of AGNs is comparable to the 
color of the standard accretion disk or slightly redder. 
The plots of all target AGNs in the (B - V) to (V - I) 
diagram show that redder targets are distributed in the 
direction of the reddening vector, starting from the bluest 
color of variable component in the diagram, which sup- 
ports Winkler's idea of the existence of an intrinsic AGN 
color. From statistical discussion of straight-line flts, it 
is found that the scatter around the best-flt regression 
line in the {V,I) plots tends to be larger than that in 
the {B, V) plots. The origin of such large scatter only in 
the {V, I) plot most likely could be attributed to the con- 
tamination from dust torus emission in the / band. All 
of these results will make rapid progress in constraining 
the mechanism of AGN variability in the future. 
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Fig. 1. — Observed light curves in the B, V, I, and J bands for the NGC3227 nucleus from 2001 November to 2007 July. 
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Fig. 2. — Observed light curves in the B and V bands for the NGC3516 nucleus from 2005 January to 2007 June. 
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Fig. 3. — Observed light curves in the B, V, I, and J bands for the NGC4051 nucleus from 2001 March to 2007 July. 
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Fig. 4. — Observed light curves in the B, V, I, and J bands for the NGC4151 nucleus from 2001 March to 2007 August. 
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Fig. 7. — Observed light curves in the B, V, I, and J bands for the MrkllO nucleus from 2003 February to 2007 August. 



8 



>;6 

X 

3 



4- 



2- 



Mrk590 


1 ' 1 




« *• ^ 


% * ' }♦ 




J-13 
'-6 ■ 




^ ^< ^ 


en V— 2 




a1 A ^'^A 
1 . 1 


A B - 



52000 



53000 



54000 



MJD 

Fig. 8. — Observed light curves in the B, V, I, and J bands for the Mrk590 nucleus from 2001 September to 2007 August. 



16 



Sakata et al. 



16 

14 

12 

?10 
E 

J 8 
6 
4 
2 



Mrk817 




* • 

* . % 


- 
■ 




1+1 


e© ^ 


« \ V+1 - 






$ 

A ^ ^ 

1 . 1 


A ~ 
A A B— 1 



53000 



54000 



MJD 



Fig. 9. — Observed light curves in the B, V, I, and J bands for the Mrk817 nucleus from 2003 September to 2007 August. 
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Fig. 11.— Observed light curves in the B and V bands for the PG0844+349 nucleus from 2003 February to 2007 May. 
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Fig. 12.— The result of surface brightness fitting for the HST images of NGC3227, NGC3516, NGC4051 and NGC4151. Shown are 
the HST/AGS images corrected for the geometric distortion {left panels), the model images of host galaxy+PSF {second left panels), the 
residual images calculated by subtracting the host-galaxy model from the nucleus-free images {second right panels), and the nucleus-free 
images calculated by subtracting the distortion-uncorrected model PSF from the original HST/ACS images then correcting them for the 
geometric distortion {right panels). All images are aligned as the north is up and the east is left. 
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Fig. 13.— The result of surface brightness fitting for the HST images of NGC4593, NGC5548, MrkllO and Mrk590. Others are same as 
in Figure 12. 
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Fig. 15. The result of surface brightness fitting for the MAGNUM V-band images of NGC3227, NGC3516, NGC4051 and NGC4151. 
Shown arc tiic original MAGNUM V-band images (left panels), tlic model images of host galaxy+PSF {second left panels), the residual 
images calculated by subtracting the host-galaxy model from the nucleus- free images {second right panels), and the nucleus- free images 
calculated by subtracting the model PSF from the original images {right panels). All images are aligned as the north is up and the east is 
left. 
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Fig. 16. The result of surface brightness fitting for the MAGNUM V-band images of NGC4593, NGC5548, MrkllO and Mrk590. 
Others are same as in Figure 15. 
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Fig. 17.— The result of surface brightness fitting for the MAGNUM V-band images of Mrk817, 3C120 and PG0844+349. Others are 
same as in Figure 15. 




Fig. 18. — The comparison between the host galaxy flux estimated from MAGNUM y-band images and that from HST F550M-band 
images, for our 11 targets. The diagonal line represents the line of agreement between the two estimates. 




Fig. 19. — The comparison between the host galaxy fluxes estimated from MAGNUM l/-band images decomposed with the galaxy-shape 
parameters being fixed at ours and with that of Bentz et al. (2009). The diagonal line represents the line of agreement between the two 
estimates. 
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Fig. 20.— The B-band flux to V-band flux diagrams for target AGNs of NGC 3227, NOG 3516, NGG 4051, NGG 4151, NGC 4593 and 
NGC 5548. Open circles represent the data, from which the HOST+NL component is not subtrax:ted, and stellate symbol represents the 
HOST+NL component. Thick line represents the best-fit linear regression. Bottom panel shows the residual plot of the data from the 
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Fig. 22.— The V-band flux to /-band flux diagrams for target AGNs of NGC 3227, NOG 4051, NGG 4151, NGC 5548, Mrk 110 and 
Mrk 590. The same as in Figure 20. 
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Fig. 23.— The same as in Figure 22, but for Mrk 817 and 3C120, 



16 

^14 
>. 

X 10 
3 



8^ 
6 
4 



2 

10 ^ 
2^ 
<]-10 |- 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I' I I I I 

3C120 




1 1 1 1 n 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ 



2 4 6 8 10 12 14 16 18 20 22 

/ flux (mJy) 



30 



Sakata et al. 




Optical Color Variability of Nearby AGNs 



18 



I I I I I I I I I I I I I I I I I I I I I I I I I I I y> 



NGC3227 




1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 



1 |- not corrected 
2^ H 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



>° 1 r corrected 
2^ ^ 
<]-10 r 



I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I 



18 20 22 24 26 28 30 32 

/ flux (mjy) 

T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

60 -NGC4151 



1 50 I- 



X 



40 - 



30 - 



Co 10 
2^ 



<1 -10 ■- 



Co 10 
2^ 



<l -10 r 




Y\^ \ I I I I I I I I I I I I I I I I I I I I I i 



'- not corrected . 



: | I I I I I I I I I I I I I I I I I I I I I I I | : 

r correctegl 



J I I I I I I I I I I I I I I L 



J I I I I I L 



30 



40 



50 



60 



70 



/ flux (mJy) 



14 - 



-12 



X 

3 



10 - 



8 - 



-I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I A , I 

NGC4051 




I l ^n I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I I I I 



vo 1 |- not corrected 
2^ H 

<]-10 r 



>° 1 r corrected 
2^ ^ 
<]-10 r . 



. 4 



_l I I I I I I I I I I I I I I I I I L. 



16 18 20 22 24 

/ flux (mJy) 



7 - 



5 6 



5 - 



Co 10 

2:^ 

<1-10 I- 



-I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 

NGC5548 




I I I I I I I I I I I I I ii I I I 



— I — I — I — I — I — I — I — I — I — I — I — I — I — |-T1 

not corrected^ ♦ * ^ 



C5 1 1- corrected 
2:^ 
<l-10 



I I I I I I I I I I I I I I I I I I I 



_l I I I I I I I I I I I I I I I I I L. 



31 



10 



11 



12 



/ flux (mJy) 



Fig. 25. — The V-band flux to /-band flux diagrams before and after the correction for dust torus emission for target AGNs of NGC 
3227, NGC 4051, NGC 4151 and NGC 5548. Open circles and dashed line show the data and the best-fit linear regression, respectively, 
before the correction. Filled squares and thick line show the data and the best-fit linear regression, respectively, after the correction. The 
stellate symbol represents the non-variable HOST-I-NL component. Middle and bottom panels show the residual plot of the data from 
the best-fit linear regression before and after the correction for dust torus emission, respectively. The residual here is expressed as a ratio 
relative to the mean observed flux for which the HOST-t-NL component has been subtra<;ted. 
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Fig. 26.— The same as in Figure 25, but for Mrk 110, Mrk 817 and 3C120. 
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Fig. 27.— The B-band flux to i?-band flux diagrams for PG0026+129, PG0052+251, PG0804+762, PG0923+201, PG0953+415, and 
PG1545+210, taken from Giveon et al. (1999). These are luminous targets in their sample of PG quasars. Thick line represents the best-fit 
linear regression. Bottom panel shows the residual plot of the data from the best-fit regression line. The residual here is expressed as a 
ratio relative to the mean observed flux. 
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Fig. 28. — The power-law spectral index of AGN continuum plotted against rest-frame frequency. Open symbols are our results and 
filled symbols are the mean colors of variable component of some Seyfert galaxies by Winkler (1997). The solid horizontal line represents 
the value of a,/ = 1/3 expected from the standard accretion disk (Shakura & Syunya«v 1973). The broken and dash-dotted horizontal lines 
represent the fitted index of QSO composite spectrum in the UV region and in the optical region, respectively (Vanden Berk et al. 2001). 
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Fig. 29. — The {B — V) to {V — I) diagram for AGN continuum. Open symbols are our results and filled symbols are the results by 
Winkler (1997). The arrow indicates the reddening vector of Ay = 1 from the bluest color of intrinsic AGN color presented in Winkler 
(1997). 
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TABLE 1 
Characteristics Of Target AGNs " 



Object 


R.A. 


Dec 


Redshift 


Galaxy size 


Av 


Morpholoj 


;y 




(J2000.0) 


(J2000.0) 




(') 


(mag) 


Host 


AGN 


NGC3227 


10h23m30.6s 


+ 19d51m54s 


0.003859 


5.4 X 3.6 


0.075 


SAB/pcc 


Syl.5 


NGC3516 


Ilh06m47.5s 


+72d34m07s 


0.008836 


1.7 X 1.3 


0.140 


SBO 


Syl.5 


NGC4051 


12h03m09.6s 


+44d31m53s 


0.002336 


5.2 X 3.9 


0.043 


Sbc 


Syl.5 


NGC4151 


12hl0m32.6s 


+39d24m21s 


0.003319 


6.3 X 4.5 


0.092 


Sab 


Syl.5 


NGC4593 


12h39m39.4s 


-05d20m39s 


0.009000 


3.9 X 2.9 


0.082 


SBb 


Syl 


NGC5548 


14hl7m59.5s 


+25d08ml2s 


0.017175 


1.4 X 1.3 


0.068 


so/a 


Syl.5 


MrkllO 


09h25ml2.9s 


+52dl7mlls 


0.035291 


0.7 X 0.4 


0.043 


So 


Syl 


Mrk590 


02hl4m33.5s 


-00d46m00s 


0.026385 


1.1 X 1.0 


0.124 


Sa 


Syl. 2 


Mrk817 


14h36m22.1s 


+58d47m39s 


0.031455 


0.6 X 0.6 


0.022 


SBc 


Syl.5 


3C120 


04h33mll.ls 


+05d21ml6s 


0.033010 


0.8 X 0.6 


0.986 


SO/BLRG 


Syl 


PG0844+349 


08h47m42.4s 


+34d45m04s 


0.064000 




0.122 


Sa 


Syl 



The data are taken from NED except for the morphology of Mrk 110 and PG0844 (Bentz et al. 2009) 
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TABLE 2 

Observational Parameters Of Target AGNs by MAGNUM telescope 



Object Filter Observed period Sampling interval^ Number of observations 











B 


V 


I 


B 


V 


/ 


NGC3227 


BVIJ 


2001.11.13 


- 2007.07.01 


11 


6 


27 


62 


123 


17 


NGC3516 


BV 


2005.01.16 


- 2007.06.26 


22 


14 




17 


35 




NGC4051 


BVIJ 


2001.03.14 


- 2007.07.30 


10 


7 


38 


67 


158 


13 


NGC4151 


BVIJ 


2001.03.14 


- 2007.08.07 


13 


5 


36 


78 


219 


22 


NGC4593 


BV 


2005.01.21 


- 2007.07.09 


20 


13 




22 


36 




NGC5548 


BVIJ 


2001.03.23 


- 2007.08.20 


10 


3 


46 


109 


302 


25 


MrkllO 


BVIJ 


2003.02.06 


- 2007.06.11 


27 


14 


134 


24 


56 


7 


Mrk590 


BVIJ 


2001.09.11 


- 2007.08.08 


8 


6 


32 


58 


128 


26 


Mrk817 


BVIJ 


2003.09.04 


- 2007.08.18 


44 


13 


122 


18 


64 


9 


3C120 


BVIJ 


2003.02.09 


- 2007.08.09 


26 


11 


34 


29 


62 


14 


PG0844+349 


BV 


2003.02.23 


- 2007.05.20 


135 


32 




9 


22 





''Median value of sampling interval (days) 
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TABLE 3 

HST/ACS/HRC Archival Images Used'' 



Object 


Exposure time'' 
(sec) 


Proposal ID"^ 




ouu 


yoo -L 


NGC3516 


300 


10516 


NGC4051 


300 


9851 


NGC4151 


300 


9851 


NGC4593 


300 


10516 


NGC5548 


600 


9851 


MrkllO 


600 


9851 


Mrk590 


600 


9851 


Mrk817 


300 


9851 


3C120 


300 


9851 


PG0844+349 


120 


9851 



''All images were obtained through the F550M 
filter. 

''Adopted exposure time of archival HST/ACS 
data 

'^Proposal ID of archival HST/ACS data 
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TABLE 4 

Parameter Values Fitted To HST/ACS/HRC Images 



dc Vaucoulcurs' Exponential 



Object 


Scale length 


Axis ratio 


Scale length 


Axis ratio 


Sky level 


Reduced 


NGC3227 


2.40 


0.58 


54.0 


0.68 


2.72 


1.774 


NGC3516 


19.1 


0.73 


1.44 


0.61 


3.89 


1.318 


NGC4051 


3.18 


0.81 


43.0 


0.74 


3.73 


1.478 


NGC4151 


10.2 


0.86 


56.6 


0.80 


2.30 


1.764 


NGC4593 


17.2 


0.58 


2.22 


0.75 


3.46 


1.347 


NGC5548 


5.65 


0.91 


6.52 


0.81 


3.48 


1.184 


MrkllO 


6.23 


0.91 






2.58 


1.143 


Mrk590 


11.6 


0.95 






5.46 


1.241 


Mrk817 


0.51 


0.67 


3.90 


0.75 


4.65 


1.119 


3C120 


0.079 


0.32 


1.31 


0.92 


5.21 


1.113 


PG0844+349 


7.7 


0.678 






0.25 


1.190 



Note. — 



Scale length in units of arcsec, and sky level in units of counts. 
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TABLE 5 

Host Galaxy Flux In The = 8". 3 Aperture 



Object 


B 




V 




/ 




Flux 


N 


Flux 


N 


Flux 


N 


NGC3227 


4.18 ± 0.32 


7 


8.60 ± 0.32 


7 


20.83 ± 0.54 


7 


NGC3516 


8.63 ± 0.19 


3 


16.07 ± 0.28 


5 






NGC4051 


4.12 ± 0.28 


5 


8.05 ± 0.38 


7 


16.84 ± 0.78 


5 


NGC4151 


9.31 ± 0.37 


7 


18.69 ± 0.90 


7 


30.05 ± 0.80 


8 


NGC4593 


3.53 ± 0.07 


4 


7.26 ± 0.094 


5 






NGC5548 


2.07 ± 0.08 


6 


4.23 ± 0.04 


6 


8.57 ± 0.14 


7 


MrkllO 


0.38 ± 0.07 


7 


0.735 ± 0.087 


6 


1.57 ± 0.11 


3 


Mrk590 


1.99 ± 0.13 


8 


4.75 ± 0.07 


9 


10.26 ± 0.09 


7 


Mrk817 


0.46 ± 0.03 


3 


1.16 ± 0.07 


8 


2.59 ± 0.29 


5 


3C120 


1.89 ± 0.32 


7 


3.50 ± 0.73 


8 


5.93 ± 1.11 


3 


PG0844+349 


0.42 ± 0.02 


2 


0.69 ± 0.16 


7 







Note. — 



Flux in units of mjy. N is the number of images used for fitting. 
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TABLE 6 
Host Galaxy Color 



Object 




B 




V 


V - I 


Hubble type 






(mai 


?) 


(mag) 




NGC3227 


0, 


.90 


± 


0.09 


1.37 ± 0.05 


SAB 


NGC3516 


0, 


,80 


± 


0.03 




SBO 


NGC4051 


0, 


.85 


± 


0.09 


1.21 ± 0.07 


Sbc 


NGC4151 


0, 


.88 


± 


0.07 


0.92 ± 0.06 


Sab 


NGC4593 


0, 


.90 


± 


0.03 




SBb 


NGC5548 


0, 


.89 


± 


0.04 


1.17 ± 0.02 


so/a 


MrkllO 


0, 


.84 


± 


0.24 


1.23 ± 0.14 


Sc 


Mrk590 


1, 


.06 


± 


0.07 


1.24 ± 0.02 


Sa 


Mrk817 


1, 


.12 


± 


0.10 


1.28 ± 0.13 


SBc 


3C120 


0, 


.79 


± 


0.30 


0.98 ± 0.27 


SO 


PG0844+349 


0, 


.65 


± 


0.26 




Sa 



Note. — Typical color of bulge is B — V = 1.02 (Kinney 
ct al. 1996), and typical color of galaxies \s (B — V,V — I) = 
(0.96, 1.31) for E, (0.78, 1.21) for Sab, and (0.50, 1.07) for 
Sod (Pukugita et al. 1995). 
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TABLE 7 
Narrow Line Flux 



Object 


[O III]A5007 


[O III]A4959 


H/3 


H7 


Reference 


NGC3227 


693 


231 


61 


28 


1, 


2 


NGC3516 


557 


186 


56 


26 


1 




NGC-1051 


408 


136 


41 


18 


1 




NGC4151 


12800 


4320 


1097 


471 


1, 


2 


NGC4593 


188 


63 


19 


9 


1 




NGC5548 


623 


208 


66 


30 


1, 


2 


MrkllO 


230 


77 


16 


7 


1, 


2 


Mrk590 


60.4 


20.1 


13.9 


6.3 


1, 


2 


Mrk817 


133 


44 


8 


4 


1, 


2 


3C120 


819 


273 


86 


42 


3 




PG0844+349 


60.3 


20.1 


6.1 


2.8 


4 





References. — (1) Whittle et al. (1992); (2) Kaspi et al. (2005); 
(3) Tadhunter et al. (1993); (4) Miller et al. (1992). 

Note. — Flux in units of 10~^^ergs cm~^. All the data tabu- 
lated here are corrected for the Galactic extinction. The [O III]A4959 
flux and H7 flux are assumed to be specific fractions of the [O III]A5007 
flux and H/3 flux, respectively. 
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TABLE 8 

Flux Contribution Of Narrow Lines In Each Filter 



Object 


Filter 


[O III]A5007 


[O III]A4959 




H7 


Total flux 


NGC3227 


B 


0.150 


0.062 


0.02 


0.019 


0.256 




V 


0.502 


0.077 


0.000 


0.000 


0.579 


NGC3516 


B 


0.100 


0.046 


0.020 


0.018 


0.184 




V 


0.474 


0.104 


0.000 


0.000 


0.577 


NGC4051 


B 


0.092 


0.038 


0.017 


0.013 


0.159 




V 


0.277 


0.035 


0.000 


0.000 


0.312 


NGC4151 


B 


2.80 


1.18 


0.44 


0.32 


4.74 




V 


9.10 


1.31 


0.00 


0.00 


10.41 


NGC4593 


B 


0.033 


0.015 


0.007 


0.006 


0.061 




V 


0.160 


0.035 


0.000 


0.000 


0.196 


NGC5548 


B 


0.081 


0.039 


0.020 


0.022 


0.162 




V 


0.606 


0.170 


0.005 


0.002 


0.783 


MrkllO 


B 


0.014 


0.007 


0.003 


0.005 


0.030 




V 


0.248 


0.079 


0.012 


0.004 


0.343 


Mrk590 


B 


0.006 


0.003 


0.004 


0.005 


0.016 




V 


0.063 


0.019 


0.006 


0.002 


0.090 


Mrk817 


B 


0.010 


0.005 


0.002 


0.003 


0.019 




V 


0.141 


0.045 


0.005 


0.002 


0.193 


3C120 


B 


0.056 


0.028 


0.019 


0.031 


0.135 




V 


0.874 


0.279 


0.059 


0.023 


1.235 


PG0844+349 


B 


0.001 


0.000 


0.000 


0.002 


0.004 




V 


0.067 


0.022 


0.007 


0.000 


0.096 



Note. — 



Flux in units of mjy. 
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TABLE 9 

Statistics Of Linear Flux To Flux Fit 



Object 


Paired bands 




6^ 


N 


Reduced 


Rv'(xlOO) 


NGC3227 


BV 


0.96 ± 0.02 


-5.26 ±0.22 


54 


0.88 


1.70 




VI 


0.78 ± 0.03 


-7.58 ± 0.82 


17 


3.30 


4.76 


NGC3516 


BV 


1.02 ± 0.02 


-8.37 ± 0.44 


9 


0.87 


1.43 


NGC4051 


BV 


0.99 ± 0.02 


-4.08 ±0.23 


55 


1.52 


1.85 




VI 


0.96 ±0.04 


-8.46 ± 0.79 


9 


1.90 


4.94 


NGC4151 


BV 


1.03 ±0.01 


-16.0 ±0.1 


36 


8.47 


0.63 




VI 


0.79 ±0.01 


3.16 ±0.25 


19 


28.53 


2.03 


NGC4593 


BV 


1.05 ±0.01 


-4.05 ±0.14 


16 


1.92 


1.67 


NGC5548 


BV 


1.01 ±0.01 


-2.60 ±0.04 


81 


1.65 


0.88 




VI 


0.70 ±0.02 


-1.14 ±0.22 


20 


5.74 


5.47 


MrkllO 


BV 


0.94 ±0.01 


-0.60 ±0.04 


11 


3.58 


0.97 




VI 


0.92 ±0.02 


-0.53 ±0.10 


7 


2.91 


4.71 


MrkSQO 


BV 


0.97 ±0.06 


-2.85 ±0.32 


19 


1.05 


5.51 




VI 


0.77 ±0.07 


-3.08 ± 0.73 


17 


4.27 


10.70 


Mrk817 


BV 


0.93 ± 0.03 


-0.66 ±0.19 


8 


0.87 


1.44 




VI 


1.07 ±0.05 


-2.53 ±0.42 


6 


0.80 


2.43 


3C120 


BV 


1.11 ±0.02 


-3.40 ± 0.26 


25 


1.70 


1.27 




VI 


1.04 ±0.01 


-2.49 ±0.23 


11 


8.05 


1.98 


PGOSll+319 


BY 


0.87 ±0.03 


-0..3S ± 0.21 


8 


1.90 


2.11 



Note. — AT is the number of data used for fitting. Ry is the ratio of fitting residual to 
flux variation defined as Rv = a-y^ /AV where cr^^ is the standard deviation of the residuals 
from the best-fit regression line, and Ay = fy"^^ — /™'° is the amplitude of flux variation. 

"fY = ax fx + b. {Xy<\x) 
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TABLE 10 

Linear V Flux to I Flux Fit With and Without Dust Correction 



With correction Without correction 



Object 


fHNa. 
J J 






avi 


Reduced 


OiVI 


Reduced 


NGC3227 


47.9 


1700 


-0.6 


-0.64 


1.52 


-0.67 


3.30 


NGC4051 


37.8 


1950 


0.0 


-0.08 


0.64 


-0.11 


1.90 


NGC4151 


51.9 


1650 


-0.6 


-0.66 


4.48 


-0.65 


28.5 


NGC5548 


16.4 


1700 


-0.2 


-0.57 


1.34 


-0.97 


5.74 


MrkllO 


1.96 


2050 


-0.6 


-0.42 


1.35 


-0.23 


2.91 


Mrk817 


5.73 


1850 


-0.7 


0.13 


0.15 


0.19 


0.80 


3C120 


7.30 


1850 


-0.3 


-0.24 


1.79 


0.10 


8.05 



''Flux of HOST+NL component in the J band in units of mjy. 
''Dust temperature in units of K. 
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TABLE 11 

Flux Contribution of the Host + Narrow-Line 

Component 



y!HOST+NL j 



Object 




B 








V 






I 




NGC3227 


0.495 


± 


0.037 


0, 


.622 


± 


0.023 


0.726 


± 


0.019 


NGC3516 


0.564 




0.012 


0, 


.714 


± 


0.014 








NGC4051 


0.482 


± 


0.032 


0, 


.643 


d= 


0.030 


0.745 


± 


0.034 


NGC4151 


0.437 


± 


0.017 


0, 


.668 


d= 


0.023 


0.649 


± 


0.017 


NGC4593 


0.472 


± 


0.010 


0, 


,681 


± 


0.009 








NGC5548 


0.565 


± 


0.021 


0, 


,789 


± 


0.025 


0.837 


± 


0.013 


MrkllO 


0.137 


± 


0.024 


0, 


,287 


± 


0.029 


0.297 


± 


0.021 


Mrk590 


0.878 


± 


0.056 


0, 


,920 


± 


0.014 


0.951 


± 


0.008 


Mrk817 


0.091 


± 


0.006 


0, 


,209 


± 


0.013 


0.322 


± 


0.036 


3C120 


0.183 


± 


0.029 


0, 


,368 


± 


0.060 


0.393 


± 


0.073 


PG0844+349 


0.077 


± 


0.004 


0, 


,118 


± 


0.024 









Note. — jHOST+NL/javc ^Yic contribution of the host- 
galaxy plus narrow-line component to the average flux of the 
light curve data. The error comes from that of the HOST-I-NL 
flux. 
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TABLE 12 

Color and Spectral Index Of AGN Continuum 



Object 


B-V 


0!BV 


V -I 


avi 




(mag) 




(mag) 





NGC3227 


0, 


.16 


± 


0, 


.02 


-0, 


.17 


± 


0, 


.07 


0, 


.66 


± 


0, 


.04 


-0, 


.64 


± 


0.11 


NGC3516 


0, 


.10 


± 


0, 


.02 


0, 


.08 




0, 


.08 




















NGC4051 


0, 


.13 


± 


0. 


.02 


-0, 


.05 


± 


0. 


11 


0, 


.37 


± 


0. 


.05 


-0. 


.08 


it 


0.12 


NGC4151 


0, 


.09 


± 


0, 


.01 


0, 


.13 


± 


0. 


,03 


0, 


.67 


± 


0. 


.01 


-0. 


.66 


± 


0.02 


NGC4593 


0, 


.06 


± 


0, 


.02 


0, 


.24 


± 


0, 


,08 




















NGC5548 


0, 


.11 


± 


0, 


.01 


0, 


.04 




0, 


,04 


0, 


.63 


± 


0. 


,04 


-0. 


,57 


± 


0.10 


MrkllO 


0, 


.19 


± 


0, 


.02 


-0, 


.32 


± 


0. 


.09 


0, 


.57 


± 


0. 


,02 


-0. 


,41 


± 


0.06 


Mrk590 


0, 


.15 


± 


0, 


.07 


-0, 


.12 


± 


0. 


.31 


0, 


.69 


± 


0. 


.20 


-0. 


.72 


± 


0.50 


Mrk817 


0, 


.19 


± 


0, 


.03 


-0, 


.32 


± 


0. 


13 


0, 


.35 


± 


0. 


.05 


0. 


.13 




0.14 


3C120 


0, 


.02 


± 


0, 


.02 


0, 


.45 




0, 


,08 


0, 


.50 


± 


0. 


.02 


-0. 


.24 


± 


0.04 


PG0844+349 


0, 


.27 


± 


0, 


.05 


-0, 


.65 


± 


0. 


,24 





















Note. — Color and spectral index are derived from equations 7 and 
8, respectively, in the text. 
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TABLE 13 
Comparison Of Internal Extinction*' 

Object Our estimate Literature Reference 

NGC3227 0.16 ±0.02 0.18 1 

NGC4151 0.090 ± 0.008 0.091 ± 0.007 2 

NGC5548 0.109 ±0.008 0.073 ± 0.018 2 

3C120 0.017 ±0.018 0.417 ±0.630 2 

References. — (1) Crenshaw et al. (2001); (2) 
Paltani et al. (1996). 

''Extinction is the value of Eb-v- 



